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ABSTRACT 

We propose a new approach for measuring the mass profile and shape of groups and clusters of 
galaxies, which uses lensing magnification of distant background galaxies. The main advantage of lensing 
magnification is that, unlike lensing shear, it relies on accurate photometric redshifts only and not galaxy 
shapes, thus enabling the study of the dark matter distribution with unresolved source galaxies. We 
present a feasibility study, using a real population of z > 2.5 Lyman Break Galaxies as source galaxies, 
and where, similar to galaxy-galaxy lensing, foreground lenses are stacked in order to increase the signal- 
to-noise. We find that there is an interesting new observational window for gravitational lensing as a 
probe of dark matter halos at high redshift, which does not require measurement of galaxy shapes. 

Subject headings: cosmology, dark matter, dark energy, gravitational lensing 
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1. INTRODUCTION 

Dark matter halos can be used as a probe of both cos- 
mological parameters and structure formation. Their sta- 
tistical distribution as a function of mass and redshift is a 
strong probe of dark energy (Allen ct al. 2004), and a com- 
parison of the halo mass versus galaxy distribution within 
individual halos contains important clues a bout the role of 
dark matter in the baryonic mass buildup ([Hoekstra et al.l 
|2005() . Any study using dark matter halos, statistically 
or individually, requires an estimate of the mass profile 
or the halo mass. Unfortunately, there are still practi- 
cal difficulties with these measurements, particularly for 
high redshift halos and certain mass range. High redshift 
clusters, for instance, are dynamically young objects, mak- 
ing velocity dispersion and X-ray temperature based mass 
estimates questionable. There is therefore a vivid inter- 
est in finding a reliable mass calibration tool: one that 
is unbiased and, as much as possible, independent of the 
dynamical state of the halo. In fact, the mass calibra- 
tion issue is central to high reds hift cluster searches, such 
as the SpARCS/SWIRE survey (iMuzzin et al.l [20091) and 
Sunyaev-Zel'dovich surveys. 

The most prevalent technique for estimating halo masses 
uses gravitational lensing, and more specifically the tan- 
gential shear measured from the shapes of background 
galaxies. The variation of the strength of the shear from 
the halo center is an estimate of the lens mass profile. Al- 
though this technique is formally unbiased and provides 
a mass estimate independent of the halo dynamical state 
and shape, it suffers from practical limitations. Accurate 
galaxy shape measurement requires good pixel resolution 
and a small Point Spread Function (PSF), which, for high 
redshift sources, is generally not the case. Current ground- 
based surveys show that it is unlikely one can reliably mea- 
sure the shape of galaxies at redshifts higher than z ~ 1.5. 
This upper bound automatically sets a maximum redshift 
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for the lenses one can probe: roughly half of the source red- 
shift. One concludes that it is unlikely shear observations 
can probe the dark matter distribution at redshifts higher 
than ~ 1. Note that, in principle, space-based data 
should perform better than ground-based data, although 
a quantitative analysis of space-based over ground-based 
performance at these redshifts remains to be done under 
realistic observing conditions. 

The alternative mass estimate we propose is based on 
lensing magnification instead of shear. Magnification re- 
lies solely on photometry instead of shapes, therefore any 
practical requirement on how large a source galaxy must 
be vanishes; technically one could measure lensing mag- 
nification on unresolved galaxies, well beyond the red- 
shift limits imposed on shape measurement by the PSF. 
Magnifica tion behind da r k matter halo s was pioneered by 
iBroadhurst et al.l ()1995[ ): iTavlor et al.l ([l9M), but at that 
time, optical surveys were too small and shallow, severely 
limiting photometric redshifts estimates. Moreover, the 
net loss in signal-to-noise with magnification compared to 
shear m ade the for mer less interesting for practical appli- 
cations (jSchneider et al. 2000 ). 

On the other hand, recent wide and deep surveys are well 
suited for measuring the magnification signal. The mea- 
sure ment of cosmic magnificat ion on Lyman Break Galax- 
ies ([Hildebrandt et al]|2009bl ) already demonstrates that 
the weak lensing signal on z > 2.5 galaxies is detectable. 
This paper explores this new opportunity for the study of 
cluster/group-sized dark matter halos at high redshift. 

Section 2 introduces the magnification technique. Sec- 
tion 3 uses the population of high redshift galaxies de- 
tected in the Canada- France-Hawaii Telescope Legacy Sur- 
vey (CFHTLS) to demonstrate how accurately halo masses 
and shapes can be measured. The last Section discusses 
the prospects for this new technique in light of current and 
forthcoming surveys. 

2. GALAXY-HALO LENSING IMAGNIFICATION 

The gravitational lensing of the source galaxies is de- 
scribed by the amplification matrix. A, which is a function 
of shear, 7, and convergence, k (see iMunshi et all (|2008[) 
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for a review). The magnification, fi{6), along a line-of- 
sight quantifies the change in flux of a distant object af- 
fected by gravitational lensing: 
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(1) 



where 6 = (0i,02) is the position angle on the sky mea- 
sured from a chosen reference point. The apparent mag- 
nitude, m, of a lensed object is m -t- 2.51og(/x), mak- 
ing it brighter (/z > 1) or fainter (/i < 1) depending 
on the amount of projected mass present in the direc- 
tion 6, as indicated in equation [1] The magnification 
effect is measured by counting, for lines-of-sight at dif- 
ferent 6, the number of galaxies that appear in a magni- 
tude bin [m, m -f dm] for a survey of limiting magnitude 
miiin > m. The number of unlensed galaxies in [to, m+dm] 
is NQ(m)dm and the number of lensed galax ies in the same 
magnitude range is N{m, 9)dm. .Naravaui (|l989t ) showed 
that: 



N{m,e)dm ^ ^2.5s(m)-i^^(^)^^^ 



(2) 



where s(rn) is the slope of the logarithmic galaxy number 
counts at magnitude m. The magnification, fj,, can be di- 
rectly obtained from a measurement of the galaxy counts 
contrast, SN{d): 
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Behind a dark matter halo, galaxies are expected to be 
magnified because both k and 7 are positive. A halo could 
also be located in an under-dense region, leading to a dim- 
ming of the background galaxies locally; however, with the 
average lensing effect over the entire sky being zero, this 
environmental effect vanishes and only contributes to the 
noise. Therefore, we can safely focus on the mass profile 
of a halo itself, embedded in a uniform background. 

The strategy adopted in this work is to evaluate 
the constraints on dark matter halos described by the 
universal Navarro-F renk- White (NFW) profile, pnfw 

avarro et al.lll997D . The mass density profile is de- 
scribed by the simple formula: 
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where pcrit is the critical density, is a characteristic scale 
radius and is the density parameter: 
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3 ln(l + C200) - C2oo/(l + C200) 
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iiW = 2 ('/'II ~ "^22) 

72(0) = 4^12 

i^{0) essentially measures the projected mass along the 
line-of-sight at 0: 
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The concentration parameter, C200, at the virial radius, 
7-200, is simply given by C200 = 7'20o/?'s- The para- 
metric model in equation 2] effectively depends on only 
two parameters: the concentration, C200j and the halo 
mass, M200, both within the virial radius and related by 
-^^200 = 4/37rr|QQPcrit- The velocity dispersion is defined 

as V^2oo = A/200 /7'2oo- 

K and 7 are given by the second order derivatives of the 
gravitational lensing potential, V'(^)- 

^ Note that the most massive clusters generate a strong magnification sij 
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where Ejvfvk(^) is the projected mass of the NFW pro- 
file, and Doi, Dqs, Dig are the angular diameter distances 
between the source^ observer and the lens. Analytical ex- 
pression of the lensing potent i al for the NFW profile can 
be found in iMeneghetti et all (|2003[ ). One could m prin- 
ciple constrain the parameters C200 and V200 for each halo 
where the magnification is measured. In practice however, 
the signal-to-noise per halo is too low and it is necessary 
to stack the signal from many foreground halos, a strat- 
egy similar to galaxy-galaxy lensing. This technique is 
described in detail in the following section. 

3. HALO CHARACTERIZATION 

Galaxy-galaxy lensing was initially developed as a probe 
of galactic-size dark matter halos, making use of the tan- 
gential shear around for eground galaxies (|Brainerd et al.l 
1 1996t [Hudson et alii 19981 ). Similar to the shear, the signal- 
to-noise ratio of the lensing magnification is low for most 
lenses, making the stacking of foreground lenses necessary 
in order to lower the noise of the average magnification 
as a function of distance from the lenses centers. The 
stacking should rely on a proxy to group foreground lenses 
with similar mass (e.g. stellar mass). It can be applied to 
clusters, groups of galaxies and individual galaxieqj- Con- 
straints on the mass and concentration parameters, M200 
and C200, of the average magnification profile are obtained 
from the likelihood: 

C a exp [{5^9) - ~SN{0))CsXi^N{0) - ^iv(^))^] (7) 
where Sn{6) is the galaxy count profile from equation 
[21 circularly averaged for our purpose, and Sn{0) is the 
galaxy count profile model. The covariance matrix, Cs^^Sk , 
is est imated as de s cribed below. 

In iTavlor et"all (|1998[ ). it was shown that for a given 
population of source galaxies, the net signal-to-noise loss 
for magnification is larger by roughly a factor of 5 com- 
pared to shear measurements. Ignoring the sampling vari- 
ance and the effect of large scale structures, there are two 
sources of noise for magnification; 1) the statistical, Pois- 
son, noise due to discrete sampling of the source galax- 
ies and 2) the clustering of the source galaxies leading 
to variations of number counts coherent over large angu- 
lar distances. In order to test the feasibility of our ap- 
proach under realistic observing conditions, the sampling 
and clustering sources of noise in Cs^Sn are measured from 
a real population of source galaxies. This population of 
source galaxies is the high redshift, z > 2.5, Lyman Break 
Galaxies (LBGs) selected f rom the C F HTLS-Wide d ata 
by the method described in iHildebrandt et al l (j2009a|) : it 
contains ^ 130 000 u- and g- dropouts in the magnitude 

jnal where stacking is not necessary l|Tavlor et al.|[T998IV 
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range 23.5 < i < 24.5 selected from 156 sq. deg. of tlie 
CFHTLS-Wide dataset. 



(CFHT i-band magnitude S 24.5; z < 2.5; clusters S 5 galaxies) 
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Fig. 1. — Cluster number counts from the Millen nium Simulation 
mock light cones of IKitzbichler and Whiti ll2007h as a function of 
redshift for different mass bins. The counts have been normalized to 
a 50 sq. deg. survey. The depth of the galaxy catalogue matches the 
CFHTLS-Wide data with i < 24.5. Clusters are defined as detailed 
in the text. 

The slope of the number counts is such that 2.5 s{rn) — 1 = 
1.4, as measured on the CFHTLS-Deep data with a lim- 
iting magnitude app roximately 2 magnitudes deeper than 
the CFHTLS-Wide (|Hildebrandt et al.ll2009bD . 

The estimate of the noise covariance matrix, Cs^Sn i pro- 
ceeds as follows: approximately 200 random positions are 
chosen to represent foreground lenses in each square de- 
gree of the 156 sq. deg. Their angular positions are then 
cross-correlated to the LBGs for each square degree. The 
average cross-correlation is zero, and the dispersion around 
zero corresponds to the Csj^Sn expected for 200 foreground 
lenses. The amplitude of the noise covariance matrix is 
later adjusted to the actual number of stacked halos. Our 
procedure takes into account realistic sampling and clus- 
tering noise of the source population. 

The realistic number of dark halos to stack for a 
given mass was take n from the mock light cones of 
IKitzbichler and White! ()2007[ ). which were created from 
the Millennium Simulatio n (ISpringel et aL 



20051) de - 



rived galaxy catalogues of iDe Lucia and BlaizotI ( 20071 ) . 
Clusters in thes e light cones are defined according to 
iMilkeraitis et al.l (|2010( ): cluster members are identified 
as having the M200 flag of their parent halo, as well as 
the same friends-of-friends identification number to avoid 
dark matter substructures. Clusters are required to have 
at least 5 galaxy members and the light cones are filtered 
out to a limiting magnitude of ?McgaCam = 24.5, to be in 
line with the CFHTLS data discussed herein. 

Figure [T] shows the number of clusters per redshift for 
a ground-based survey of 50 sq. deg.: for redshifts > 1 
we expect ~ 10^ clusters of mass ~ IO^^Mq and ~ 150 
clusters of mass > IO^^A/q. Note that the power spec- 
trum normalisation in the Millennium Simulation (MS) is 
fJs = 0.9 and thu s substantially large r than the WMAP7 
value, as = 0.8 (|Larson et al.ll2010l) : therefore, the MS 
over-predicts the number of halos by a factor of approxi- 
mately 3 — 4. The cluster number counts shown in Figure 
[l]are thus approximately the expected number in our Uni- 
verse for a survey of ~ 200 sq. deg. with the CFHTLS- 
Wide depth. 

The performance of the magnification measurement is 



evaluated for a low and a high mass bin, using halo counts 
taken from Figure [TJ The low mass bin consists of 6000 
stacked halos with ^200 — 650 km/s and C200 = 5-5 and 
the high mass sample contains 135 stacked halos with 
V200 = 950 km/s and C200 = 4.5. All halos are placed 
at redshift z = I and we assume no centroid misalignment 
in the stacking. Figure [5] shows the constraints on V200 
and C200 for mass bins using the magnification measured 
on a source of LBGs as defined previously, assuming they 
are all at redshift z — 3. For a fixed concentration, the 
velocity dispersion for both mass bins is constrained with 
an accuracy better than 5%. This shows that magnifica- 
tion of high redshift galaxies can in principle probe dark 
halos well beyond the domain of applicability of the shear 
method. By extrapolating, it is clear that precision mea- 
surement of dark halo mass profiles is possible with future 
full sky lensing surveys. 

The practical implementation of this technique will re- 
quire a mass proxy and an estimate of the cluster/group 
center which coincides with the halo peak. The mass 
proxy is not specified at this stage, but it has to rely 
on indicators such as the stellar mass. X-ray temper- 
ature/luminosity, cluster richness, or a combination of 
these. I nterestingly, the dispersio n of the mass-richness 
relation (iHilbert and White! !2010! ) appears significantly 
larger than the velocity dispersion shown in Figure [2j This 
strongly suggests that the calibration of halo richness as a 
mass proxy using lensing magnification can be refined for 
even smaller stacked foreground samples. A centroid mis- 
alignment would cause a spread in the mass profile, as sug- 
gested in !Mandelbaum et a"l~l (j2006l) and !Sheldon et al.! 
(200'3) . A complete study of the impact of center misiden- 
tification is beyond the scope of this paper. Note that 
for the two extremes of halo sizes, individual galaxies and 
massive clusters, the centroiding issue is nearly nonexis- 
tent. 

An interesting aspect of magnification is the possibil- 
ity of simultaneously measuring the lensing signal and the 
dust absorption by the lenses. Dust absorption by the lens 
induces a small chromatic angular scale-dependent cross- 
correlation between the background population and the 
lens (unlike le nsing magnification which is achromatic). 
!Menard et al.! (2009a, b) showed that the dust absorption 
is a very negligible source of magnitude noise and that 
its contamination is only a few percent of the lensing sig- 
nal. Only small angular scales are affected, where bary- 
onic matter is concentrated and dust absorption is highest. 
The practical implementation of how dust absorption can 
be fully integrated into a magnification based mass mea- 
surement is left for a forthcoming study. 

4. CONCLUSION 

We have shown that weak lensing magnification of z > 
2.5 Lyman Break Galaxies can reach a relatively high 
signal-to-noise to enable the study of dark matter halos 
at high redshift, z > 1, where traditional shear measure- 
ments would fail because the source galaxies are unre- 
solved. The performance of our method was quantified 
using a real LBGs distribution, from the CFHTLS-Wide 
data, as source galaxies. The steep number counts and 
the high redshift of the sources considerably help alleviat- 
ing the low signal-to-noise found in previous magnification 
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Fig. 2. — Top panels: the left plot shows the velocity dispersion, V20O1 versus concentration parameter, c: constraints from the stacking 
of 135 halos with V200 = 950 km/s and C200 = 4.5 at redshift 2 = 1; on the right, the constraints from the stacking of 6000 halos with 
V200 = 650 km/s and C200 = 5.5 are shown. Bottom panels: left and right show the halo mass derived from the constraints (V20O1 c) from the 
top left and top right panels respectively. 



studies. The approach is similar to galaxy-galaxy lensing, 
where lenses are stacked according to a mass proxy, al- 
though no attempt was made to discuss which mass proxy 
should be used since it entirely depends on the data wave- 
length coverage. 

The magnification technique as a probe of dark matter 
halos can be generalized to a larger sample of background 
galaxies over a wider magnitude and redshift range. It 
is interesting to note that, in principle, the halo shape 
can also be measured, which can be used to discriminate 
between various Cold Dark Matter models. The main lim- 
itations of our method lie in our ability to separate lensing 
from dust extinction in large surveys, which is dependent 
on the number and wavelength coverage of filters, and in 
identifying reliable mass proxies and determining the halo 
centers, although the latter two equally affect shear based 



mass profile measurements, it is not specific to magnifica- 
tion. A quantitative estimate of the practical limitations 
for future lensing surveys such as LSST, and JDEM would 
be particularly interesting, but is beyond the scope of this 
work. 
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